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We calculate the density of states and the the thermal conductivity k in layered superconduc- 
' tors for rotating in-plane magnetic field by using approximate analytical expressions for Abrikosov's 
I vortex lattice state. For a gap with d^^y - orbital symmetry we find that the differences AA'^ and Ak 
between field directions parallel to the antinodes and parallel to the nodes of the gap are positive 
^ I for lower fields and change sign for higher fields near Hc2 ■ For increasing impurity scattering, AA'' 
, and Ak decrease due to a broadening of certain peaks in the angle-resolved density of states. The 
^ ■ frequency dependencies of the amplitudes AA'' and Ak exhibit similar crossover behavior at low 
' frequencies which give rise to a sign change of Ak at low T/Tc ■ We conclude that variations of 
the specific heat and thermal conductivity can be observed only in the clean and low temperature 
limits. 
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C/3 _ Recently fourfold oscillations of the specific heat have been observed for rotating in-plane magnetic field in pure 
Sr2Ru04 at very low temperatures.^ These data have been analyzed on the basis of the Bogoliubov-de Gennes (BdG) 
equations and the Pesch-approximation for the vortex state with the result that the spin-triplet superconductor 
Sr2Ru04 has a gap with nodes in the directions of the a and b axes.^ In an earlier measurement of the thermal 
' conductivity k no appreciable angular variation was observed from which it was concluded that no vertical line nodes 
^ , exist in Sr2Ru04 However, the amplitude of the oscillation of k for a state with vertical gap nodes was found to be 
Q ' strongly temperature dependent and to change sign for increasing temperature.^ This explains the failure to observe 
O ] an angular variation of k because the measuring temperature in the experiment of Ref. 3 was too high. 

In the present paper we repeat the calculations of the density of states N and the thermal conductivity k for 
rotating in-plane field for a superconductor with d^^j, symmetry of the orbital part of the gap using the analytical 
J> . expressions of Pesch^'^ based on the quasiclassical Green's functions which are easier to handle than the original 
' expressions based on the Gorkov equations and Abrikosov's vortex state.''''* Our main purpose is to investigate the 
, effect of impurity scattering. This depends sensitively on the phase shift for impurity scattering as has been shown 
■ in a recent paper for d-wave pairing in fields parallel to the c-axis.^ We shall only present results for the unitary limit 
of impurity scattering with phase shift 7r/2 because the relevance of the unitary limit has been strongly suggested by 



the results for the universal heat transport in Sr2Ru04 ^ 
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^ For simplicity we neglect the small c-axis dispersion of the cylindrical Fermi surface and thus the c-axis component 

. of the Fermi velocity. The azimuthal angle of the quasiparticle momentum in the ab-plane with respect to the a-axis 

■ is denoted by (f> and and we use a for the direction of the magnetic field in the ab-plane. Thus the component of the 

, Fermi velocity in the direction perpendicular to the magnetic field is given by v sin((/) — a) . For the momentum part 

I of the order parameter we assume d^y symmetry which, together with the field dependence of the spatial average of 

' O the modulus of the Abrikosov vortex state, yields the expression 

3: |Ap = A^[l-;^2]V2sin2(20) (1) 

K*" ' where h = H/Hc2 and H is the spatial average of the field. The Pesch-approximation now yields for the spatial 
' average of the density of states N{uj) , normalized to the normal state density of states A'o : 

C3 ; N{Lj)/No = Re{l + 8\A\^[A/vsm{<j)~a)]^[l + iy/^zw{z)]y^^^ ; (2) 

where 

z = 2[w + iSj][A/t;|sin(0-a)|] ; A^{2eH)-^/^. (3) 

Here is the self energy for impurity scattering which is calculated self-consistently in the t-matrix approximation 
from the expressions given in Ref. 8. The thermal conductivity k in the BPT-approximation is given in Ref. 4 and, in 
the equivalent Pesch-approximation in Ref. 6. The integrand of the w-integral in Ref. 6 is proportional to the density 
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of states [Eq.(2)] times the lifetime of the quasiparticles. This hfetime is equal to the reciprocal of the sum of the 
scattering rates due to impurity and Andrecv scattering which have been discussed at length in Ref. 8. 

In Fig.l we show our results for N{lu = 0)/Na and k(T = 0)/k„ (k„ is the normal state thermal conductivity) versus 
the reduced field h = H/Hc2 for fields along the gap node (a = , solid curves) and along the antinode (a = 7r/4 , 
dashed curves). The heat current for k is along the antinode of the gap, (p = t^I^ ■ The impurity scattering self energy 
is calculated self-consistcntly for the unitary phase shift limit and the reduced scattering rate 5 = F/Ap = 0.01 where 
r = l/2Tn is the normal state scattering rate. This is the right order of magnitude for the purest samples of Sr2Ru04 
with the highest One sees that the dashed curve for N/No (antinode) lies above the solid curve (node) for h from 
zero to about the crossing field h,. = 0.63 while, above he , it lies below the solid curve. The difference AA^ between 
the two curves corresponds to the amplitude of the measured fourfold anisotropy of the specific heat.^ This amplitude 
is about 6 percent of N{a = 0) at /i = 0.3 . Our results agree approximately with the results obtained in Ref. 2 by 
solving the BdG equations and by using the Pesch-approximation. For n/nn (lower curves in Fig.l) the dashed curve 
(antinodal direction of the field) lies above the solid curve (nodal direction of the field) for h between and about 
he = 0.95 where it crosses the solid curve. As with the density of states, the difference Ak between the dashed and 
solid curve corresponds the amplitude of the fourfold anisotropy and amounts to about 25% of K{a = 0) . It has been 
shown in Ref. 4 that the amplitude Ak of the fourfold oscillation in rotating field decreases rapidly with increasing 
temperature and reverses its sign at a very low temperature indicating that the variations with field direction can 
only be observed at very low temperatures. 

In addition to the temperature effect, impurity scattering also diminishes the amplitudes of the fourfold oscillations 
of the density of states and the thermal conductivity. An example is shown in Fig. 2 where we present our results for 
the unitary phase shift limit and a much larger scattering rate 6 = 0.1 which is of the order of magnitude corresponding 
to lower Tc 's in Sr2Ru04 Comparison of Fig. 2 with Fig.l shows that, for the larger scattering rate, the crossing 
field he ~ 0.4 for the curves of N is much smaller and the relative amplitude AN/N{a = 0) ^ 2% is also much smaller. 
The crossing field for n and the relative amplitude An/nla = 0) are of the same orders of magnitude as those for N . 

In Ref. 2 the reasons for the crossover behavior of the density of states have been explained in detail by calculating 
the angle-resolved density of states for low and high fields using the BdG equations and the Pesch-approximation. 
Similar angle-resolved densities of states N{(j))/No are shown in Fig. 3 for a low field h = 0.2 and nodal field direction 
(a = , solid curves) and antinodal field direction (a = 7r/4 , dashed curves). Note that, in the latter case, the angle (p 
is measured from the field direction a = -k/A. One sees that, for a = , a single broad peak occurs at <^ = 7r/2 while, 
for a = 7r/4 , two peaks occur at i/i = 7r/4 and (j) = 37r/4 . This explains why in this field region the total density of 
states N{a = it/A) is larger than N{a = 0) because, for a = 7r/4 , all four nodes contribute to N while, for a = , two 
of the four nodes are parallel to the field and are thus unable to contribute to A^ . In the high magnetic field region, 
on the other hand, N{a) = 0) > N{a = 7r/4) because then the behavior is determined by N{(l)) for ^ = and ^ = tt . 
Comparison of Fig. 3(a), calculated for the unitary limit and 5 = 0.01 , with Fig. 3(b), calculated for S = 0.1 , shows 
that the main effect of stronger impurity scattering is to broaden the peaks of the solid curve at = and (p = n and 
to enhance the minima at (p — tt/A and 37r/4 . At the same time the minima of the dashed curve (p — and = tt are 
enhanced. This effect of impurity scattering is similar to the effect of finite c-axis dispersion of the cylindrical Fermi 
surface which was taken into account in Ref. 2. The total effect of the broadening of the peaks and the enhancement 
of the minima in the angle-resolved density of states for low fields is the decrease of the amplitudes of variation of N 
and K in rotating in-plane field as shown in Figs.l and 2. 

Finally we consider the frequency dependence of N{uj)/No for nodal and antinodal field directions and its effect on 
the frequency and temperature dependence of the thermal conductivity. In Fig. 4(a) we have plotted our results for 
N{lo)/Nci versus = w/Aq for a low field {h = 0.2) for nodal (a = , solid curve) and antinodal (a = 7r/4, dashed 
curve) field directions. In Fig. 4(b) we show the corresponding curves for K{u))/Kn which is the factor multiplying 
(w/T)^ sech^(a;/2T) in the normalized integral over d{Lo/T) for This function of O = w/Aq yields approximately 

the temperature dependence of through the relation O ~ 2.4(T/Ao) because the integrand is strongly peaked 

at uj/T = 2.4 . The expression for is given by N{w)t{uj) where the quasiparticle scattering rate l/2T(ti') is equal 
to the sum of the impurity and Andreev scattering rates. ^ The most interesting results in Figs. 4(a) and 4(b) are the 
crossovers of the dashed and solid curves occuring at very low frequencies which correspond to a sign change of the 
amplitudes of the variations of fourfold symmetry with rotation angle. For k,{lj) this first crossover occurs at about 
ilc — 0.1 which means that, at a temperature of about T/Tc ~ 0.1,, the amplitude of the variation with rotation 
angle vanishes or changes sign. It should be pointed out that for high fields the curves for N{(jj) and k(w) look quite 
different. The solid curve for N{ijj) starts out at a much higher value than the dashed curve and the coherence peak 
for the solid curve vanishes. This is similar to the results of Ref. 2. 

In sumary, we have calculated the density of states A'' and the thermal conductivity k for a gap with Axy - orbital 
symmetry and magnetic fields in the directions H \ \ gap node (a = , solid curves) and H \ \ gap antinode (a = 7r/4 , 
dashed curves). For both calculations we have used the approximate expressions of Pesch^'^ for Abrikosov's vortex 
lattice state which are easier to handle than the original expressions of the BPT-approximation.^'^ Comparison of 
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these results with the rcsuhs of the quasiclassical Eilenbergcr equations^° and the BdG-cquations^ has shown that 
these analytical expressions provide very good approximations over the whole field range from i?c2 down to H^i ■ We 
have included impurity scattering and used a phase shift of 7r/2 because the unitary limit in Sr2Ru04 is strongly 
suggested by universal heat transport.^ 

Our main result is that the density of states N{ui = 0) and the thermal conductivity k{T = 0)//v„ are somewhat 
larger for the antinodal field direction (a = it /A) than for the nodal field direction {a = 0) in a field region < h < he 
{h = H/Hc2) where he denotes the field where the dashed curves cross the solid curves (see Fig.l). Our results for 
N agree essentially with the results obtained from the BdG equations and the Pesch-approximation.^ These results 
have led to the conclusion that the observed minima and maxima of the specific heat for rotating in-plane field in 
Sr2Ru04 are due to vertical line nodes in the directions of the a- and b-axes.^ It turns out that the amplitudes 
AA'' = A''(a = 7r/4) — N{a = 0) and Ak = K{a = 7r/4) — K,{a = 0) and the crossover field he decrease for increasing 
impurity scattering (see Fig. 2). The reason can be seen from a comparison of the angle-resolved density of states 
N{(f>)/NQ in Figs. 3(a) and 3(b). The most important contributions to the density of states [(see Eqs.(2) and (3)] 
are the terms A/v\ sm{(j) — a)| where A = (2ei/)^^/^ is the magnetic length, v \ sm{(j) — a)| is the component of the 
Fermi velocity perpendicular to H , and 4> and a are the azimuthal angles of the quasiparticle momentum and field 
direction, respectively. Asymptotic expansion of the w- function'' in Eq.(2) shows that in the limit cj) ^ a the density 
of states takes on the BCS form for |A(0)|^ [see Eq.(l)]. For quasiparticles moving in the directions perpendicular 
to the field {(j) — a = ±7r/2) the analytic expression involving thr w-function yields the full effect of the superfluid 
flow and the Andreev scattering due to the complex order parameter of Abrikosov's vortex lattice function. These 
directional terms, together with the angular dependence AQsin^(2(/)) of the gap function, yield quite different angle- 
resolved densities of states N{(j))/NQ for a = and a = 7r/4 shown in Fig. 3. Comparison of Figs. 3(a) and 3(b) shows 
that the main effect of stronger impurity scattering is to broaden the peaks of the solid curve at <^ = and </> = tt . It 
is interesting that this effect is quite similar to the effect of a finite c-axis dispersion of the cylindrical Fermi surface 
which has been taken into account in the analytic Pesch-expression in Rcf. 2. Finally, we have investigated the finite 
temperature effect on the amplitude Ak for the variation with rotating in-plane field by calculating N{lu)/No and 
K{uj)/K,n for the field directions a = and a = -k/A (see Figs. 4(a) and 4(b) for low field h = 0.2 and small impurity 
scattering 5 = 0.01 ). The large differences between N{u) and k{uj) arise from the quasiparticle lifetime t{(jj) in the 
latter expression which is given by the reciprocal of the sum of scattering rates due to impurity and Andreev scattering. 
The impurity scattering rate is calculated self-consistently in the unitary limit of the t-matrix approximation, and the 
Andreev scattering rate arises from the imaginary part of the self energy proportional to A(ri)A*(r2)G(ri — r2, — w) , 
where A(r) is Abrikosov's vortex lattice order parameter^ and G is the hole propagator. The crossover of the dashed 
and solid curves at a low frequency w ~ O.IAq in Fig. 4(b) leads to a corresponding crossover in the temperature 
dependencies of k, at about T/Tc — 0.1 showing that the amplitude of the variation with field rotation vanishes and 
changes sign at this temperatiirc. Similar results have already been obtained in Ref. 4 using the BPT-approximation. 

In conclusion we can say that the amplitudes of the variations of the specific heat and thermal conductivity in 
rotating in-plane field for a gap with vertical line nodes decrease rapidly with increasing impurity scattering and 
increasing temperature. 
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FIG. 1. 1) Density of states, N{u! = 0)/No , and thermal conductivity, K/Kn{T — » 0) , versus h = H/Hc2 for in-plane field 
direction parallel to the gap node (a = , solid curves) and parallel to the antinode (a = 7r/4, dashed curves). The reduced 
impurity scattering rate is 5 = F/Ao = 0.01 with phase shift in the unitary limit. 
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Fig. 2 




Fig. 3a 




FIG. 3. 3a) Angle-resolved density of states, N{4>)/No for h = 0.2. Solid curves are for nodal field direction {a — 0) and 
dashed curves for antinodal direction (a = 7r/4). The angle 4> is measured from the field direction. Impurity scattering with 
5 = 0.01 in the unitary limit. 



7 



Fig. 3b 




FIG. 4. 3b) Angle-resolved density of states, N{4>)/No for h = 0.2. Solid curves are for nodal field direction (q = 0) and 
dashed curves for antinodal direction (a = 7r/4). The angle 4> is measured from the field direction. Impurity scattering with 
5 = 0.1 in the unitary limit. 
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Fig. 4a 




FIG. 5. 4a) Density of states N{uj)/No versus fl = lo/Aq for h = 0.2 , 5 = 0.01 , and field directions a = (solid curve) and 
a = 7r/4 (dashed curve). 
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Fig. 4b 




FIG. 6. 4b) Thermal conductivity K{u!)/Kn versus f2 = uj/Aq for h = 0.2 , S = 0.01 , and field directions a = (solid curve) 
and a = 7r/4 (dashed curve). The temperature gradient is in the direction of the antinode. 
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